APS/123-QED 



The Conversion of nonlocal one-body operators into local ones: 

The Slater potential revisited 



James P. Finley 
Department of Physical Sciences, Eastern New Mexico University, 

Station #33, Portales, NM 8813^ and 
Department of Applied Chemistry, Graduate School of Engineering, 
The University of Tokyo, Tokyo, Japan 113-8656 
(Dated: February 2, 2008) 

Abstract 

One-particle Schrodinger equations are considered, e.g., the Hartree-Fock equations, that contain 
a nonlocal operator, e.g., the Hartree-Fock exchange operator, where this operator depends on the 
one-particle density-matrix of a determinantal state. One-body nonlocal operators of this type are 
converted into approximate local potentials that depend on the kernel of the nonlocal operator 
and, also, the one-particle density matrix that, as mentioned above, the nonlocal operator also 
depends on. When the non-local operator is the exchange operator, the method yields the Slater 
potential. 



I. INTRODUCTION 



There is a considerable interest in the conversion of nonlocal, one-body operators into 
local, multiplicative operators, or potentials. For example, using a local, multiplicative 
operator to replace a nonlocal one can lead to significant computational savings when solving 
the Hartree-Fock equations, or other similar one-body, coupled equations. In addition, the 
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Kohn-Sham approach of density functional theory (DFT) 



|7J requires the 



noninteracting state to be obtained from a local potential, since this formalism invokes 
the Hohenberg-Kohn theorem [^, even though, according to the adiabatic connection, the 
exact exchange energv is the same one-particle density-matrix functional as the one from 
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except that the one-particle functions are Kohn- 



Hartree-Fock theory 
Sham orbitals. 

Other possible applications of nonlocal to local conversions is for one-particle density- 
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For example, 



matrix approaches that employ a nonlocal operator 
in the reference-state one-particle density-matrix formalism 2y, |2jJ , the conversion of the 
nonlocal exchange-correlation potential into a local operator leads naturally to a density 
functional theory, since the one-particle density matrix of the determinantal reference state 
obtained from a local potential is a unique functional of its electron density, according to 
the Hohenberg-Kohn theorem. 

A classical example of the conversion of a nonlocal operators into a local one is Slater's 
local exchange-operator, derived from a uniform electr on g as (3, 13] • A very similar ap- 



proximate functional has also been derived by Gasper [2^, where this one only differs by 
a constant from the one derived by Slater. The Xa approach 2^ includes a semiempirical 
parameter that yields Slater's original prescription with a — 1, and the one by Gasper for 
a = |. Gasper's exchange operator is also used within the local density approximation 
(LDA) of DFT When the Gasper potential is combined with Becke's exchange 

correction [3], derived in 1988, a very accurate approximation is obtained, and this two 



component exchange potential 
known as the B3LYP method 



is used within the very successful DFT hybrid approach 
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The optimized potential method (OPM) 
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is an approach to 



convert a nonlocal operator into a local potential. This method, unfortunately, leads to 
rather complicated functionals that depend on the individual orbitals. Below we introduce 
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an approach that converts one-body nonlocal operators into approximate local potentials 
that depend on the kernel of the nonlocal operator and, also, the one-particle density matrix 
that the nonlocal operator also depends on. When the non-local operator is the exchange 
operator, the method yields the Slater potential j^, l^jj which also appears as the 
leading term from the Krieger-Li-Iafrate (KLI) approximation of the OPM [29 
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II. THE CLOSED-SHELL RESTRICTED SPIN-ORBITAL FORMALISM 

Consider a set of spatially restricted spin-orbitals [38]: 

ip ia (r h u) = x;(ri)crOi), a = a, /3, (1) 

where the spin and spatial coordinates are given by r x and uii, respectively. Consider also 
the following spinless, one-particle Schrodinger equation: 

F pi Xi(?x) = e»X«(ri), (2) 
where this Hermitian one-body operator is given by 

K =-^V 2 + t; Pl +u; pl , (3) 

and v pi and w Pl are local and nonlocal operators, respectively; these operators depend on 
the spinless one-particle density matrix pi of a closed shell single-determinantal state, say 
|$); that is, we have [J HQ 



Pi 



> 1 ,r 2 ) = 2^ Xw (r 1 )x;N 1 (4) 



where the same spatial-orbitals appearing in this summation are doubly occupied within 
|$); henceforth, w, x denote these occupied orbitals from |$); r, s denote the unoccupied 
orbitals. 

Now consider the possibility of replacing the nonlocal operator w pi by a local one, say z pi ; 
so, we have 

z Pl (ri)xx{ri) = J dv 2 u; pi (ri,r 2 )x*(r2), (5) 
where w Pl (ri, r 2 ) is the kernel of w Pl . 
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Our interest here is in the one-particle density- matrix p\ that arises from solving Eq. (J2J) . 
Therefore, we only consider the operator w Pl acting upon an occupied orbital, w Pl x x in the 
above equation. In other words, pi does not depend on w pi x r , where Xr is an excited orbital. 
We do, however, require w pi to be Hermitian. 

Multiplying the previous equation by Xx( r z) anc ^ summing over the orbital indices gives 

z Pl {Y 1 )pi{xi,v z ) = j rfr 2 w Pl (ri,r 2 )pi(r 2 ,r 3 ). (6) 

Setting r 3 = ri yields the desired result: 

z Pi( r i) ~ P~ l { T l) J dr 2 w pi (r 1 ,r 2 )p 1 (r 2 ,r 1 ), (7) 

where p(r) is the electron density, pi(r, r). For example, the exact exchange operator from 
Hartree-Fock theory yields the Slater potential 22, 3, 36], 

<(ri) « —p-'in) J dv 2 r u 1 \p 1 (r 1 ,v 2 )\ 2 , (8) 

where the kernel of the exchange operator is — |rj~ 2 1 pi(ri, r 2 ). 

In the two previous expressions above we have changed the equality to an approximation, 
since these are apparently not identities, where this conclusion arises, in part, since, in our 
derivation of Eq. © we have summed over all occupied orbitals. Now if Eq. (jSj) is valid, 
of course, Eqs. © and (J7j) must be satisfied, but, not vice versa. So if we define z pi by 
Eq. ((7j) we will probably not satisfy Eq. (j3J). For example, consider a simple case of only 
two occupied orbitals, say Xw and Xx, where instead of Eq. (jSJ) being satisfied, we have the 
following relations: 

^(rOXzOi) = / rfr 2 w Pl (ri,r 2 )xx(r 2 ) + x^(ri)0(ri), (9a) 
z P i( r i)Xy(ri) = J dv 2 w pi (r 1 ,r 2 )x y (r 2 ) - ^(r 1 )0(r 1 ), (9b) 

where 0(r 3 ) is an arbitrary function. Multiplying the first equation by Xx( r 3)* and the 
second one by Xy( r 3)* removes the last terms from both equations when the two equations 
are added, yielding Eq. (P). 

As in any operator, our operators are completely defined by their matrix elements. If the 
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following identities are satisfied, then the local and nonlocal operators are equivalent: 

J dr i Xy{ r i)z P i( r i)Xx{ri) = J dri J dv 2 x y (ri)w Pl (ri, r 2 )Xx(r 2 ), (10a) 
J dri Xr(ri)* pi (ri)xs(ri) = J dr t J dr 2 Xr(r 1 )w pi (r 1 ,r 2 )xx(r 2 ), (10b) 
J dri Xr(ri)« w (ri)x*(ri) = J dr x J dr 2 Xr{ri)w pi (r h r 2 )x s (r 2 ), (10c) 

and we can replace the one-body operator, given by Eq. (jSJ), by the following: 

^i = ~V 2 + « w +^. (11) 

However, for our purposes, we do not need all three relations given by Eqs. (fTUJ) to be 
satisfied. In particular, if Eq. ()10c|) is not satisfied, Eq. (J2J) is still satisfied for the same 
occupied orbitals and, as mentioned previously, the one-particle density-matrix pi is not 
changed. Furthermore, if Eq. (jlUaj) is not satisfied, we will get different occupied orbitals 
but they will differ only by a unitary transformation as long as Eq. (llObJ) remains valid, and, 
again, the one-particle density-matrix pi is not changed. Hence, we only need Eq. (jlObjl to 
be a reasonable approximation. 

III. THE OPEN-SHELL UNRESTRICTED SPIN-ORBITAL FORMALISM 

We now generalize the previous derivation to the case where the determinantal state is 
composed of orbitals that are spatially unrestricted; that is, we have 

^(xi) = Xia(ri)a(uj), a = a,P, (12) 

where the spatial and spin coordinates U\ are denoted collectively by x x ; furthermore, we 
have a different one-body operator for each spin-function: 

F°Xia(ri) = eiaXia(ri), (13) 



f; = --v 2 + ^ + ^;, (i4) 



where the Hermitian one-body operators are given by 

1, 

2 

and v° and w° are local and nonlocal operators, respectively; these operators depend on the 
one-particle density matrix 7 of the single-determinantal state |$); that is, we have 

7(x 1; x 2 ) =7(ri,u;i;x2,Wi)<$ W iua, (15) 



where 

7(ri,w;x 2 ,u;) = ^ \ cr ( UJ )\ 2 Pia(ri,r 2 ), (16) 
and for iV CT occupied a orbital, the spin-components of 7 are given by 

Pi<x(ri, r 2 ) = ^ Xwa(ri)x* wa (r 2 )- (17) 

Our one-body operators F°, and potential, and w° depend upon 7. However, 7 is an 
explicit functional of its components, p\ a and p\$\ furthermore, we require F° not to depend 
on the spin variable u, since this dependence is easily removed. So, we could, instead write, 
for example, F° la Pl/3 ■ However, for a less cluttered notation, we will continue to indicate a 
7 dependence. This is also not inaccurate, since if we know 7, we also know p\ a and p\p. 

Now, as in Eq. (jSJ), consider the possibility of replacing the nonlocal operator w° by a 
local one, say z": 

z°(ri)Xxa(ri) = J dr 2 w°(r 1 ,r 2 )Xx*(r2), (18) 

where w°(ri, r 2 ) is the kernel of w°. Multiplying this equation by x* xa { v ^) an d summing over 
the orbital indices gives 

2^(ri)pi CT (ri,r 3 ) = J dr 2 w°(ti, r 2 )pi CT (r 2 , r 3 ). (19) 

Setting r 3 = ri yields the desired result: 

4( ri ) w p- 1 ^) J dr 2 ^(r l5 r 2 )p 1(T (r 2 , r x ), (20) 

where p CT (r) is the cr-component of the electron density, given by pi CT (r, r). For example, the 
exact exchange operator from Hartree-Fock theory yields the Slater potential |22l. I23L l36||. 

t^(n) « -p; 1 ^) I dr 2 r^l^Cn,^)! 8 , (21) 

where the kernel of the exchange operator is — rj" 2 1 p lcr (r 1 , r 2 ). 
APPENDIX A: THE SLATER POTENTIAL 

The Slater potential is a local approximation of the exchange potential that was suggested 
by Slater and used to obtain the Xa approach with (a = 1). While it is well known in the 
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literature that Eq. (|21|) is the Slater potential, we are not familiar with a derivation demon- 
strating that is the case. Therefore, we have decided to present one here. In appendix IB1 we 
present a derivation of the Xa exchange potential. 



Using the notation by Slater 
form: 



22| , but in atomic units, Slater's potential has the following 



Vs Xi = " ' n T( 7 , r , Al 

where the sums are over the N occupied spin orbitals from the set {uj}. Using our spin- 
orbitals, given by Eq. (JT2J), we have 

, ) = E^7^( x i)^(x 2 )^(xi)^ CT (x 2 )rr 2 1 rfx 2 

E.Ef ^(xi)^(xi) 

where the sums are over the N a and Np occupied spin orbitals, and from Eq. (fT2*J) . we have 
u.(xi) = (A3) 

Y. N x °Y.„xlM)xUri)\°M\ 2 

Using the following identity: 

ffViWl) = Saa'^Ml 2 , (A4) 

gives 



v s (xi) = ^ — - , (A5) 

E<t l°"(^i)l 2 E 



which can be rearranged, 



«.W = — 7 „ ' V \ 1 , (A6) 

E>K)l 2 (Ef ^(ri)x* CT (rx) 



E. 




Pi a (r 1 ,r 2 )\ 2 r 1 ^ dr 2 




o-K)IV(ri) 



and by using Eq. (fT7|) . we have 

..<*> - -£i;a;.-- (at) 

Furthermore, the following identity is readily verified: 

f s (ri,o;i)xi CT (ri)cr(a;i) = tj* CT (ri)xi CT (ri)cr(u;i) (A8) 

indicating that the Slater potential v s , given by Eq. (jAl|) or ()A2|) . is equivalent to our 
potential, , given by Eq. (}2"T]) . 
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APPENDIX B: THE Xa EXCHANGE POTENTIAL 



It is well known that the spinless, one-particle density matrix pi(ri,r 2 ) of a closed-shell 

unified electron-gas depends on its density p, a constant. A generalized expression for pi is 

n 

obtained by replacing the constant density with a non-constant one, say p(r) |3[: 

/ \ / \ „ / \ /sin \kf(r)s] — kf(r)scos \kf(r)s]\ , 
Pi(ri, r 2 ) = Pl (r, s) = 3p(r) I {J^M 3 ~^^J ' (B1) 

where 



k f (r) = ^3vr2p(r), (B2) 

and a change of coordinates is employed: 

r = ^(r 1 + r 2 ), (B3) 
s = ri -r 2 , (B4) 
s = |s| = r 12 = [(x 2 - zO 2 + (3/ 2 - yi ) 2 + (z 2 - Zl f] 1/2 . (B5) 

Using this expression for p 1; the exchange-energy functional for a closed shell system, given 
by 

E x [pi] = ~J J dridr 2 r^ 1 \p 1 (v 1 ,r 2 )\ 2 , (B6) 
leads to the well known Dirac exchange-energy density-functional 

Ex[p] "~l(-) 1/3 / p4/3(r)dr ' (B7) 
where in the derivation by Parr and Yang |3j , the integration is expressed using the r and s 
coordinates and the integration is carried out over s. Analogous expressions for unrestricted 
orbitals and open shell systems are also readily derived 0, |y| . 

In the Kohn-Sham LDA ODD , the functional derivative of the above exchange func- 
tional gives the following local exchange potential: 

3 X 1/3 

-P(r) , (B8) 

7T / 



in agreement with Gasper's potential |2J]. However, it well known that this functional differs 
from the one obtained by Slater j^] by a factor of ~, where v^(r) is obtained directly from 
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the Hartree-Fock, nonlocal, exchange-operator and, in addition, Slater's approach uses a 
sort of averaging over the occupied orbital states. 

We now demonstrate that the Xa exchange potential, with original prescription of a = 1, 
arises in our approach when considering a uniform electron gas, where we only consider the 
closed-shell spin restricted formalism, and the derivation in analogous to the one by Parr 
and Yang mentioned above. In our derivation, we also use Eq. (|B1|) except that, at least for 
the moment, we leave the density ClS db constant: 

pi(r 1; r 2 ) = 3pp^ (sin(k f s) - k f s cos{k f s)) . (B9) 

Substituting this expression into Eq. (jSJ) and using Cartesian coordinates gives 

v* 1 (r 1 )=~pj J J dx 2 dy 2 dz 2 -J^g (sm(k f s) - k f scos(k f s)) 2 . (BIO) 

Now as far as the integral is concerned, x±, yi, are Z\ are constants. Therefore, by making 
the following substitution: 

s x = (x 2 - xi) (Bll) 
as well as analogous ones for s y and s z , where, for example (ds x = dx 2 ), we get 

v* pi (r x ) = ~pj J J ds x ds y ds z -^ (sm(k f s) - k f s cos(k f s)) 2 , (B12) 
which in spherical coordinates, can be written as 



-18nj- j ds —— g (sin(kfs) — kfS cos(kfs)) 2 . (B13) 
kf J s kj 



Making the following substituting t = kfS, where ds = dt/kf, we have 



Vp^rt) = -187r^2 / dt ^ (ant -t coat) . (B14) 
kf Jo t 

and since the integral is | p, we get 

where we have introduced a different notation for the exchange potential, since it is simply a 
constant that depends on p. Using Eq. ()B2jl . we obtain the following expression after some 
algebra and by replacing p by a non-constant density, p(r): 



3 /3 



1/3 



t£(r) = --(-p(r)) . (B16) 



This expression is in agreement with Slater's original prescription of a — 1. 
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